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Abstract 

Although subgrid-scale models of similarity type are insufficiently dissipative for 
practical applications to large-eddy simulation, in recently published a priori analyses, 
they perform remarkably well in the sense of correlating highly against exact residual 
stresses. Here, Taylor-series expansions of residual stress are exploited to explain the 
observed behavior and “success” of similarity models. Until very recently, little atten- 
tion has been given to issues related to the convergence of such expansions. Here, we 
re-express the convergence criterion of Vasilyev et al. [J. Comput. Phys., 146 (1998)] 
in terms of the transfer function and the wavenumber cutoff of the grid filter. 
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1 Introduction 


In contrast to direct numerical simulation (DNS), in which alt energetic scales of motion are 
resolved on a fine grid, in large-eddy simulation (LES), the spatially filtered Navier-Stokes 
equations (FNSE) are solved numerically on a relatively coarse grid. For incompressible flow, 
and in tensor notation, the FNSE are given by 


du k 

dx k 

dU k d(V k Vt) 
dt dxi 


0 

ap 1 ffgfc thu 
dxk Redxidxi dxi 


( 1 ) 

( 2 ) 


where Uk is the velocity vector, p is the pressure, and r« is the subgrid-scale 


stress 


tensor (or residual-stress tensor) defined exactly as 

7jy 5 SfcSj — ttfctti (3) 


Here overlines denote grid-filtered quantities, and Re is the Reynolds number. 

The residual-stress tensor incorporates the effects of the unresolved scales of motion upon 
the resolved scales. It is customary in LES to model these effects, for which there exists a 
variety of possible models. In 1991, Germano et al. [1] introduced the concept of dynamic 
modeling, which exploits the resolved- turbulent-stress tensor £y > a computable quantity that 
is extracted by applying an explicit secondary filter-the test filter-to the resolved velocity 
fields as follows: 

At = Afri — EfcCt (4) 

Here hats denote test-filtered quantities. In general, the test- and grid-filter widths, A and 
A respectively, may differ. For later convenience, we denote their ratio by r = A/A. 

Recently, interest in SGS models of similarity type (e.g., Liu et al. [2], Stolz and Adams 
[3], and Pruett and Adams [4]) has revived (despite the fact that these models are typically 
insufficiently dissipative for stand-alone applications to LES). Under certain conditions, re- 
markably high correlations between r k i and At have been observed in a priori analyses based 
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on experimental (Liu et al. [2]) and computational (Pruett and Adams [4]) data. On the 
basis of their observations, Liu et al. [2] propose the stress-similarity model 

Tki & ci,Cki (5) 

where ci is simply a constant. 

i •; f "E,[ ' ; ■ i : 

In this paper, we exploit Taylor-series expansions of the grid- and test-filter operators to 
analyze the relationship between the tensors t*j and Cm and to optimize the value of ci in 
Eq. 5. The convergence properties of such expansions are subtle, and attention is devoted 
herein to establish a criterion for convergence. For simplicity, we apply filtering only in 
the x direction and suppress the y and z coordinates and the time t. However, numerical 
experiments (Pruett and Adams [4]) suggest that the results carry over to multi-dimensional 
filters. 

2 Taylor-Series Analyses 

Although the use of fixed-width filters is common in LES, in our judgment, this is an ill- 
advised practice that usually leads to the contamination of the SGS dissipation by the 
truncation error of the numerical scheme. Recently, the utility of tunable (one-parameter) 
filters has been recognized (e.g., Vasilyev et al. [6] and Pruett and Adams [4]). One- 
parameter filters permit the filter width A (or preferably for this work, the wavenumber 
cutoff k c ) to be specified independently of the grid increment Ax. In our view, k c should 

' : ' ft j ■ i ' ; ; : ! ' 

be specified on the basis of physical considerations; that is, the cutoff should lie in the 
inertial range of the Kolmogorov spectrum. On the other hand, Ax should be determined 
by numerical considerations; that is, by the grid resolution necessary to resolve the smallest 
eddies for the numerical scheme of choice. The dimensionless product a c = k c Ax defines the 
cutoff parameter to remove the degree of arbitrariness. To this end, we exploit one-parameter 
filters of Pade type, the details of which are relegated to the Appendix. 

The following discussion applies to a priori analyses, in which DNS data are filtered to 
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extract the exact residual stresses and their modeled counterparts. By implication, in this 
section Ax = Axqns- As discussed in the Appendix, a discrete filter can be represented by 
a Taylor-series expansion in the grid increment Ax. For example, if u k is a filtered velocity 
field, then 

tifc(x) = Ufc(x) 4- aiti' fc (x)Ax + a 2 t4'(x)Ax 2 + 03 ttjj['(x)Ax 3 + ... (6) 

Here, primes denote (partial) derivatives with respect to x. For notations! simplicity, the 
factorials associated with the Taylor expansion have been absorbed into the coefficients a ; -. 
As shown in the Appendix, the Taylor series of a filter implies its transfer function and vice 
versa. In general, a filter is said to be of order m if the first non-vanishing coefficient of its 
Taylor series is Om. By applying Eq. 6 to Eq. 3, we obtain a Taylor series for the residua! 
stress, namely 

Tki = (of — 2o 2 )u*u{Ax 2 + (ai «2 — 3o3)(u*tt" + u|u*)Ax 3 + ... (7) 

Because the coefficients aj depend implicitly upon a c , it is useful to regard t*j — tu(q c ). As 
the tensor tm arises solely from the quadratic nonlinearity of the Navier-Stokes equations, it 
is quadratic at leading order in Ax, provided that the filter is of either first- or second-order. 
On the other hand, if the filter is of order m > 2, then is of leading order rh. 

With the help of Mathematical a similar expansion of Cm leads to 

C ki = (a 2 — 2a 2 )t4tx{r 2 Ax 2 + (ofr 2 - 2aia 2 r 2 + aia 2 r 3 — 3azr 3 )(u' k u" + ti|u£)Ax 3 + ... (8) 

where the test and grid filters differ only in their respective widths, whose ratio is r. In (pre- 
ferred) terms of grid- and test-filter cutoffs, (A: c )g and (fc c )x, respectively, r = (/? c )g/(^ c)t- 
Clearly, Cm = £fci(a c ,r). By comparing Eq. 7 and Eq. 8, we conclude that the SGS stresses 
are approximated to leading order by 



whereby cl = 1/r 2 in Eq. 5. How good is the approximation? From Eqs. 7 and 8, we obtain 
the approximation error 

Em(<x c , r) s r w - - [3a 3 (r - 1) + 0102(3 * r) - of] (v! k u" + u{«2) Ax 3 + ... (10) 
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from which we conclude the following: 


1. If the filter is of either first- or second-order in Ax, then the approximation error 
(Eq. 10) is of higher order (0(Ax 3 )) than is the subgrid-scale stress (0(Ax 2 )), and 
the approximation could be accurate provided that there is separation between the 
contributions at various orders (a topic discussed shortly). 

2. If the filter is of order two (a! = 0, 02 ^ 0) then 

Em = 3ct3(r — l)(u*u" + u{u£)Ax 3 + ... (11) 

3. Although r = 1 is precluded in LES for reasons to be addressed shortly, r = 1 is 
optimal for a priori analyses because, for second-order filters, the leading-order error 
vanishes. 

4. The approximation that results from the optimality condition r = 1 is simply the 
scale-similarity model proposed by Bardina et al. [5], namely t*/ « Uk% — ft*®/. 

Whenever the grid and test filters differ in their Taylor coefficients, the situation is 
somewhat more complicated. Here, we presume that the grid and test filters are each of 
second order and symmetric, in which case 

Uk(x) = tt*(x) + a 2 u'k(x)Ax 2 + a 4 u^(x)Ax 4 + ... (12) 

Uk(x) = u*;(x) + b2u'k(x)Ax 2 + 6 4 i4 4 *(x) Air 4 + ... 

Applying Eq. 12 to Eq. 3, we obtain 

tm = — 2 a 2 u' k u' l Ax 2 + [(a^ — 6a 4 )(uJ'u") — 4a 4 (u*u| 3 ^ + «{u^)]Aa; 4 + 0(Ax 6 ) (13) 

Similarly, from Eqs. 12 and 4, we derive 

Cki = — 2b2u'ku\Ax 2 + [( b\ - 66 4 )(t4'u") — (46 4 + 20262) (it^itj 3 * + ujitj. 3 ^) Ax 4 + 0(Ax 6 ) (14) 

A comparison of Eqs. 14 and 13 reveals the following approximation to be correct to leading 
order: 

Tkl ~ (15) 
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If the grid and test filters are identical ( bj = aj) then < 12/62 = 1, in which case most of the 
leading-order error drops out to yield 

Eki = r kl -Cu = 2a\{u' k uf ) + tiJu* 3 ) )Ax 4 + 0(Ax 6 ) (16) 

By either approach, we conclude that the use of identical grid and test filters should produce 
optimal results in a priori analyses. 

The present analysis corroborates and provides an explanation for the observations of 
Liu et al. [2] and Pruett and Adams [ 4 ] , whose a priori analyses were based on experimental 
and computational data, respectively. To be specific, Liu et al. [ 2 ] observe negligible cor- 
relations between the residual stresses and their Smagorinsky-modeled counterparts when 
using spectral filters, which act with exponential order. Our results suggest that, because the 
Smagorinsky model is of second-order in Ax, the model is appropriate only in the context of 
first- or second-order grid filters. Moreover, for spatial top-hat filters (which are of second 
order), Liu et al. [ 2 ] observe the highest correlations between r k i and C k i for r = 1 , which we 
have shown to be optimal. In a priori analyses from DNS of isotropic turbulence, Pruett and 
Adams [4] observe correlations of the form C(r k i , £«) of nearly unity whenever the filter is 
weakly to moderately dissipative (tt/2 < a c < it) and r = 1 . As the dissipation increases, the 
correlation coefficient diminishes gradually, but it remains surprisingly high (C > 0.8) even 
for quite dissipative filters (a e « 0.45). Furthermore, they observe correlations to degrade 
somewhat as r deviates substantially from unity. 

3 Taylor-Series Convergence 

Although Taylor-series analysis is frequently exploited to develop or analyze SGS-stress mod- 
els (e.g., Rogallo and Moin [7] and Horiuti [ 8 ]), the approach has sometimes been criticized 
because of uncertain convergence, an issue that has received little attention until the recent 
work of Vasilyev et al. [ 6 ]. In practice, Taylor-series approximation is most useful if the 
leading-order error is relatively small, which in turn requires that the series converges fairly 
rapidly. Here, for completeness, we draw on the work of Vasilyev et al. [ 6 ] and re-express 
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their convergence criterion in terms of the transfer function and wavenumber cutoff of the 
grid filter. 

Let a denote one of the three velocity components, and, for simplicity, presume that a 
is periodic on [— it, 7r] and contains no wavenumbers higher than fcmax> a finite integer such 
that fcmax > k v , the Kolmogorov wavenumber. It will be useful in later discussion to think 
of kmax as the wavenumber cutoff for full resolution of quadratic nonlinearities in DNS, in 
which case, for a spectral numerical scheme 

fcmaxAxoNS — n (17) 


Equation 17 is simply the Nyquist criterion of signal processing. Given the assumptions 
above, the velocity field and its derivatives may be expanded as finite Fourier series; that is, 

fcm&x fcmax 

u(x) = £ U k e ikx and u (m) (z) = JT (tjfc)" 1 ^ 1 ** f| (18) 

fcss-fcmax fc=— fcmax 

where i = >/— 1> and U k is the k-th (complex) Fourier coefficient of a. Slightly adapting 
Vasilyev et al. [6] (who also assume finite fcmax), we obtain a bound on the magnitude of 
the m-th derivative of a, namely 


fcmax 

|>* w (*)l < E l*mi 

—fcmax 



< s/iE (j*™ k^dk) 1/2 

— -J^F ( ^max \ ’ _ jLn* / Ekmax 

- - tmax v ^ttt 2 

where Holder’s inequality is used to proceed from the first to the second step of Eq. 19, and 
where E = \Uk\ 2 is the total “energy” as per Vasilyev et al. [6]. Recall from Eq. 6 


that 

u(z) = £ (oo = 1) (20) 

m=0 
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By Eq. 19 


£ \a m uW(x)Ax"'\ < £ l«m|J£x^ 5 (WAx)"* (21) 

m=0 m=0 V 771 -h 1/Z 

Thus, the series in Eq. 20 is absolutely convergent provided that the series on the right-hand 
side above (Eq. 21) converges. By the ratio test, the power series in fcmaxAx above converges 
provided that 

lim fcmaxAs <1 (22) 

m-»oo n “““• ' 


As discussed in the Appendix, the coefficients a m depend on derivatives of the transfer 
function H(a, oc c ) of the filter as follows: 


flm(Qfe) = 


(m!)i n 


From Eqs. 17, 22, and 23, we obtain the following convergence criterion in terms of the filter’s 
transfer function and its cutoff: 

«■ , ", Ig-j 1 < 1 (24) 

m-»oo ( m + 1) |Jf( m )(0,Qf c )| 

For filters whose stencils and coefficients are symmetric, the appropriate criteria are 

lim < 1 or lim T < i (25) 

By the Cauchy product theorem, if the convergence of the Taylor series for u* is guaranteed 
by Eq. 24 or Eq. 25, then the series for Uk% also converges. Moreover, by definition, &max 
is sufficiently large so that quadratic nonlinearities are well-resolved in DNS. It follows that 
if the series for ti* converges, so must that for These two additional considerations 
guarantee that the Taylor series for (Eq. 3) converges provided Eq. 24 or Eq. 25 is 
satisfied. For symmetric, fully explicit filters (see the terminology in the Appendix), Eq. 25 
holds for all values of the cutoff parameter such that the off-center weights remain of the 
same sign. 

We now specialize the analysis to the symmetric, second-order Fade filter, whose transfer 
function H(a, a c ) is given in Eq. 34 and Fig. 1 of the Appendix. For this class of filters, the 
dependence upon a e of the convergence criterion is subtle. Originally, we conjectured that 
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m 

0 

1 

2 

3 

4 

5 

0>2m 

i 

0.75 

0.4375 

0.2521 

0.1453 

0.0838 

|^2m+2|/| a 2m| 

NA 

0.75 

0.5833 


0.5766 

0.5766 


Table 1: Taylor-series coefficients for Pade filter with a c = 7r/3. 


Eq. 25 holds at least for all tt/2 < a c < n. However, attempts at formal proof failed. Subse- 
quent numerical experimentation over a range of cutoffs suggests (but does not guarantee) 
that Eq. 25 is satisfied approximately for —0.1 < C = — cos(a c ) < 1, and, as a general rule 
of thumb, the more dissipative the filter (the smaller a c ), the slower the convergence. (This 
is not entirely accurate, but a truer statement is too detailed for the present forum.) 

Here, we will be content to consider three specific cutoff values, two for which the conver- 
gence criterion Eq. 25 is satisfied, and one for which it is not. First, for DNS, o c = 7r. From 
Eq. 34, H(a,n) = 1, independently of a. Consequently, u k = u k and T k i(n) = u k U[ - u k ui = 
itjfcUi — iifcttj = 0. Similarly, for r = 1 , £*/ (tt, 1) = H/Mi - = u*uj - u k ui = u k ui - u k ui = 0. 

As expected, in the DNS limit, the filters turn off, the residual stress vanishes, the resolved 
turbulent stress also vanishes (for r = 1), and the resulting Bardina model is trivially exact. 
Second, for the special case a c = 'njl, |jf( 2m )(0, 7 t/ 2)| = 1/2 for all m > 0, in which case 
Eq. 25 is clearly satisfied. Third, we consider a c = 7r/3 (C = -0.5), for which Table 1 
presents the first few coefficients ai m and their ratios. Although the coefficients eventually 
diminish at an apparently constant rate, that rate is too slow to satisfy Eq. 25. For suffi- 
ciently dissipative Pade filters, the Taylor coefficients apparently grow at an asymptotically 
constant rate. 


Returning briefly to the (convergent) case for which ot c = tt/ 2, we have 

tm(*/ 2) = -itiiuJAi 2 - + J 2 + u i“i 3) )] Al< + - 


(26) 


and 

E kl (tt/2, 1) = r kl - £ kl = | (u' fc uj 3) + uj«i 3) ) Ax* + ... (27) 

Careful examination of Eqs. 26 and 27 reveals that the leading order error term above 
could be large relative to the magnitude of the residual stress. Indeed, criterion Eq. 25 
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guarantees only that the terms of Taylor expansions eventually diminish, not necessarily 
that the leading order approximation error is small. One might expect that higher order 
corrections to the expansion for Cm would be necessary before the similarity model (Eq. 5) 
is even approximately valid. (See, for example, the generalized similarity model of Stolz 
and Adams [3].) In this light, that correlations C(tm, Cm) are typically observed to be 
high (Pruett and Adams [4]) even for quite dissipative filters (whose Taylor series may even 
diverge) is at first intriguing. With regard to LES, there are three ameliorating factors. First, 
in the inertial subrange, the Fourier coefficients f/* decay as AT 5 / 6 in amplitude. Second, for 
fcmax > k > kr,, the Fourier coefficients decay even more rapidly due to the effectiveness 
of viscosity. Neither factor has been taken into account. Third and dominant, the bound 
expressed by Eq. 19 is extremely pessimistic; equality holds only when all Fourier components 
align in phase, an unlikely scenario in a turbulent flow. 

4 LES 


Finally, we turn to an implication of our analysis to LES. For a priori analyses, there is no 
prohibition on the use of identical test mid grid filters. However, r = 1 is disallowed in LES 
with secondary filtering, as the following line of reasoning suggests (and as also noted in Liu 
et al. [2]). For LES, Ax = Axles- Recall that (fc c ) G and (Ai c )t denote the grid- and test- 
filter cutoffs, respectively, where we now presume (k e ) t < ( k c ) G - By design, once ( k c ) G is 
established from physical considerations, then (for a spectral numerical scheme) the Nyquist 
criterion (A: C ) G Axles = (<*c)g « vr establishes the appropriate grid increment Axles- By 
definition, (& c )t Axles = (<*c)t- Thus, 


t*les = 


_ (*c)g 


> 1 


(28) 


( k c ) T (<*c)t 

In LES, it is common to use j*les = 2. Attempts to use tles = 1 yield (c* c )t = tt, which 
turns off the secondary filter. However, at least in theory, there is no reason why similarity 
models with r > 1 are not viable for LES. 
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5 Conclusions 


In conclusion, Taylor-series expansions of residual stress are lisetiil for developing and analyz- 
ing SGS models for LES. In particular, Taylor-series analysis reveals the explicit dependence 
of the residual stress on the grid filter, and thereby reveals potential model-filter inconsis- 
tencies. Moreover, such analysis underscores the desirability of one-parameter filters, whose 
cutoffs can be tuned independently of grid resolution. The degree of usefulness of Taylor- 
series expansions, however, depends ultimately on how rapidly the series converge. A rigorous 
convergence criterion has been developed that is expressed in terms of the transfer function 
and wavenumber cutoff of the grid filter. As a rule of thumb, the more (less) dissipative the 
filter, the slower (faster) the convergence. Taylor-series expansions have then been applied 
to the analysis of the similarity models of Bardina [5] and of Liu et al. [2]. The analysis 
corroborates and appears to explain a number of observations regarding such models (e.g., 
Liu et al. [2] and Pruett and Adams [4]). In a priori analyses, similarity models typically 
perform far better than expected based on the mathematical analysis herein, which suggests 
that a refined analysis is both desirable and possible. 
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Appendix: One-Parameter Pade Filters 


Following Lele [9], a discrete, second-order, one-parameter (C) Pade filter is constructed by 
considering the symmetric pointwise scheme 

C/i-i + fi + C/m — °/» + 2 (/‘ _1 **■ /*+ l ) » (* = 2> — n ~ 1) (29) 

Various treatments are possible for the boundary nodes i — 0 and $ = n. However, the 
simplest tack is to impose no filtering at the boundaries. In matrix form Aff = Ef, where 
M and E are tridiagonal. 


The action of the filter above on the single Fourier mode exp(iArx) yields the complex 
transfer function 


H(a) 


a + frcos(of) 


(30) 


1 + 2£ cos(a) 

where a = A; Ax. For applications to LES, we consider only low-pass filters, for which 
H( 0) = 1 and H(n) — 0 . In combination, these constraints imply b = a = 0.5 + C- 
Admissible values of the parameter are -5 < C < Whenever C = 0 , M is the identity 
matrix, and the filter is fully explicit. The fully explicit case corresponds to a discrete top-hat 
filter with weights at adjacent nodes of [1/4, 1/2, 1/4]. Otherwise, the filter is implicit. The 
value C = 2 yields M = E, which turns off the filter. For all admissible parameter values, 
the matrices M and E are diagonally dominant with positive diagonal elements, in which 
case the filter operator L = M~ l E is positive semidefinite. 


By analogy to discrete differentiation operators, to each discrete filter is associated a 
Taylor-series expansion of the form 

f{xi) = f(xi) + axf'(xi)Ax + a 2 /"(xj)Ax 2 4- a 3 / w (xi)Ax 3 + ... (31) 


Here, for notational simplicity, the factorials in the Taylor series have been absorbed into 
the coefficients a ; . By applying the Taylor-series representation of the filter to exp(tAxj), we 
obtain the corresponding Taylor series of the transfer function, namely 

H(a) = 1 + ai (ta) + a 2 («x) 2 + a 3 (ta) 3 + ... (32) 
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( 33 ) 


Evaluating successive derivatives of Eq. 32 at a = 0 yields 

_ HW( 0 ) 

Thus, the Taylor series implies the transfer function and vice versa. In general, a filter is of 
order m if its first non-vanishing Taylor coefficient is a*,. Filters associated with symmetric 
stencils are of even order with purely real transfer functions. In particular, the Pade scheme 
above is of second order (provided C # 0.5). 


The coefficients of the Taylor series of Eq. 31 are functions of the parameter. To quan- 
tify this dependence for the Pade filter, we (unconventionally) define a C) the dimensionless 
wavenumber cutoff, such that \ = H(a c ), whereby 


H(ot,a c ) - 


(1 - coso c )(l + cosa) 


2(1 — cos a c cos a) 

Figure 1 compares the transfer functions of the Pade filter for selected values of a c . 


(34) 
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Figure 1: Transfer function of one-parameter family of second-order low-pass filters of Pade 
type for selected values of dimensionless wavenumber cutoff a c . 
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